REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services.  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway  Suite  1204.  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law.  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  ( DD-MM-YYYY )  2.  REPORT  TYPE 

30  Sep  2009  REPRINT 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

INCORPORATING  LOVE-  AND  RAYLEIGH-WAVE  MAGNITUDES,  UNEQUAL 
EARTHQUAKE  ANDEXPLOSION  VARIANCE  ASSUMPTIONS,  AND 
INTRASTATION  COMPLEXITYFOR  IMPROVED  EVENT  SCREENING 

5a.  CONTRACT  NUMBER 

FA87 1 8-09-C-00 1 2 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6260 IF 

6.  AUTHOR(S) 

Jessie  L.  Bonner',  Anastasia  Stroujkova1,  Robert  H.  Shumway2, Robert  B.  Herrmann', 
and  Dale  N.  Anderson4 

5d.  PROJECT  NUMBER 

1010 

5e.  TASK  NUMBER 

SM 

5f.  WORK  UNIT  NUMBER 

Al 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Weston  Geophysical  Corporation 

181  Bedford  St.,  Suite  1 

Lexington,  MA  02420 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 

29  Randolph  Road 

Hanscom  AFB,  MA  01731-3010 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL/RVBYE 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AFRL-R  V-H  A-TR-2009- 1 086 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited. 

Weston  Geophysical  Corporation1,  University  of  California  Davis2,  St.  Louis  University3,  Los  Alamos  National  Laboratory4 

13.  SUPPLEMENTARY  NOTES 

Reprinted  from  Proceedings  of  the  2009  Monitoring  Research  Review  -  Ground-Based  Nuclear  Explosion  Monitoring  Technologies,  21  - 
23  September  2009,  Tucson,  AZ,  Volume  III  pp  445  -  454. 

14.  ABSTRACT 

Our  ohjective  is  to  improve  seismic  event  .screening  using  the  properties  of  surface  waves.  We  are  accomplishing  this  through  (1 )  the  development  of  a  Love-wave  magnitude  formula  that  is 
complementary  to  the  Russell  (2006)  formula  for  Rayleigh  waves  and  (2)  quantifying  differences  in  complexities  and  magnitude  variances  for  earthquake  and  explosion- generated  surface  waves 

We  have  applied  the  Ms  (VMAX)  analysis  (Bonner et  al.,  2006)  using  both  Love  and  Rayleigh  waves  to  events  in  the  Middle  Last  and  Korean  Peninsula.  For  the  Middle  Last  dataset,  consisting  of 
approximately  100  events,  the  Love  Ms  (VMAX)  is  greater  than  the  Ray  leigh  Ms  (VM  AX)  estimated  for  individual  statioas  for  the  majority  of  the  events  and  azimuths,  w  ith  the  exception  of  the 
measurements  for  the  smaller  events  from  European  stations  to  the  northeast.  It  is  unclear  whether  these  smaller  events  suffer  from  magnitude  bias  for  the  Love  waves  or  whether  the  paths,  which 
include  the  Caspian  and  Mediterranean,  have  variahle  attenuation  for  Love  and  Rayleigh  waves.  For  the  Korean  Peninsula  we  have  estimated  Rayleigh-  and  Love-wave  magnitudes  for  .11 
earthquakes  and  two  nuclear  explosions,  including  the  25  May  2009  event  For  25  of  the  earthquakes,  the  netw  ork-averaged  Love-w  ave  magnitude  is  larger  than  the  Rayleigh-wave  estimate.  For  the 
2009  nuclear  explosion,  the  Love -wave  Ms  (VMAX)  was  3. 1  while  the  Rayleigh-wave  magnitude  was  3.6. 

We  are  also  utilizing  the  potential  of  observed  variances  in  Ms  estimates  that  differ  significantly  in  earthquake  and  explosion  populations.  We  have  considered  tw  o  possible  methods  for 
incorporating  unequal  variances  into  the  discrimination  problem  and  compared  the  performance  of  various  approaches  on  a  population  of  73  Western  United  States  (WUS)  earthquakes  and  131 

Nevada  Test  Site  (NTS)  explosions.  The  approach  proposes  replacing  the  Ms  component  hy  Ms  +  a*o,  where  a  denotes  the  interstation  standard  deviation  obtained  from  the  stations  in  the  sample 
that  produced  the  Ms  value.  We  replace  the  usual  linear  discriminant  a*M.s  +h *ntb  with  a *Ms  +b *mb  +c*c  In  the  second  approach,  we  estimate  the  optimum  hybrid  linear-quadratic  discriminant 
function  resulting  from  the  unequal  variance  assumption.  We  observed  slight  improvement  for  the  discriminant  functions  resulting  from  the  theoretical  interpretations  of  the  unequal  variance 
function 

We  have  also  studied  the  complexity  of  the  “magnitude  spectra’’  at  each  station.  Our  hypothesis  is  that  explosion  spectra  should  have  fewer  focal  mechanism-produced  complexities  in  the 
magnitude  spectra  than  earthquakes.  We  have  developed  an  intrastation  "complexity"  metric,  AMs.  where  AMs  -  MsU)-Ms(i+\)  at  periods  i,  which  are  between  9  and  25  seconds.  The  complexity  by 
itself  has  discrimination  power  but  does  not  add  substantially  to  the  conditional  hybrid  discriminant  that  incorporates  the  differing  spreads  of  the  earthquake  and  explosion  standard  deviations. 

15.  SUBJECT  TERMS 

Surface  wave  magnitude.  Love  waves.  Seismic  discrimination 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

SAR 

18.  NUMBER 
OF  PAGES 

10 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Robert  J.  Raistrick 

a.  REPORT 

UNCLAS 

b.  ABSTRACT 

UNCLAS 

c.  THIS  PAGE 

UNCLAS 

19b.  TELEPHONE  NUMBER  (include  area 
code ) 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


DTIC  COPY 


2009  Monitoring  Research  Review:  Ground-Based  Nuclear  Explosion  Monitoring  Technologies 


INCORPORATING  LOVE-  AND  RAYLEIGH-WAVE  MAGNITUDES,  UNEQUAL  EARTHQUAKE  AND 
EXPLOSION  VARIANCE  ASSUMPTIONS,  AND  INTRASTATION  COMPLEXITY 
FOR  IMPROVED  EVENT  SCREENING 

Jessie  L.  Bonner1,  Anastasia  Stroujkova1,  Robert  H.  Shumway2, 

Robert  B.  Herrmann',  and  Dale  N.  Anderson4 

Weston  Geophysical  Corporation1,  University  of  California  Davis2, 

St  Louis  University3,  Los  Alamos  National  Laboratory4 

Sponsored  by  the  Air  Foree  Research  Laboratory 

Contract  No.  FA8718-09-C-0012 
Proposal  No.  BAA09-75 


ABSTRACT 

Our  objective  is  to  improve  seismic  event  screening  using  the  properties  of  surfaee  waves.  We  are  accomplishing 
this  through  (1)  the  development  of  a  Love-wave  magnitude  formula  that  is  complementary  to  the  Russell  (2006) 
formula  for  Rayleigh  waves  and  (2)  quantifying  differences  in  complexities  and  magnitude  variances  for  earthquake 
and  explosion-generated  surfaee  waves. 

We  have  applied  the  Ms  (VMAX)  analysis  (Bonner  et  aL,  2006)  using  both  Love  and  Rayleigh  waves  to  events  in 
the  Middle  East  and  Korean  Peninsula.  For  the  Middle  East  dataset,  consisting  of  approximately  100  events,  the 
Love  A/v  (VMAX)  is  greater  than  the  Rayleigh  Ma  (VMAX)  estimated  for  individual  stations  for  the  majority  of  the 
events  and  azimuths,  with  the  exception  of  the  measurements  for  the  smaller  events  from  European  stations  to  the 
northeast.  It  is  unclear  whether  these  smaller  events  suffer  from  magnitude  bias  for  the  Love  waves  or  whether  the 
paths,  whieh  inelude  the  Caspian  and  Mediterranean,  have  variable  attenuation  for  Love  and  Rayleigh  waves. 

For  the  Korean  Peninsula,  we  have  estimated  Rayleigh-  and  Love-wave  magnitudes  for  3  1  earthquakes  and  two 
nuelear  explosions,  including  the  25  May  2009  event.  For  25  of  the  earthquakes,  the  network-averaged  Love-wave 
magnitude  is  larger  than  the  Rayleigh-wave  estimate.  For  the  2009  nuelear  explosion,  the  Love-wave  Ms  (VMAX) 
was  3.1  while  the  Rayleigh-wave  magnitude  was  3.6. 

We  are  also  utilizing  the  potential  of  observed  variances  in  Ms  estimates  that  differ  significantly  in  earthquake  and 
explosion  populations.  We  have  considered  two  possible  methods  for  incorporating  unequal  variances  into  the 
discrimination  problem  and  compared  the  performance  of  various  approaches  on  a  population  of  73  Western  United 
States  (WUS)  earthquakes  and  131  Nevada  Test  Site  (NTS)  explosions.  The  approach  proposes  replacing  the  Ms 
component  by  Ms  +  a*a,  where  a  denotes  the  interstation  standard  deviation  obtained  from  the  stations  in  the  sample 
that  produced  the  Ms  value.  We  replace  the  usual  linear  discriminant  a *MS  +b *mh  with  a*A/v  +b *mh  +e*a.  In  the 
seeond  approach,  we  estimate  the  optimum  hybrid  linear-quadratie  discriminant  function  resulting  from  the  unequal 
variance  assumption.  We  observed  slight  improvement  for  the  discriminant  functions  resulting  from  the  theoretical 
interpretations  of  the  unequal  variance  function. 

We  have  also  studied  the  complexity  of  the  “magnitude  spectra”  at  eaeh  station.  Our  hypothesis  is  that  explosion 
spectra  should  have  fewer  foeal  meehanism-produeed  complexities  in  the  magnitude  spectra  than  earthquakes.  We 
have  developed  an  intrastation  "complexity"  metric,  A\fs>  where  A \fs  =  A/v(/)-A/v(/+ 1)  at  periods  /,  whieh  are  between 
9  and  25  seeonds.  The  complexity  by  itself  has  discrimination  power  but  does  not  add  substantially  to  the 
conditional  hybrid  discriminant  that  incorporates  the  differing  spreads  of  the  earthquake  and  explosion  standard 
deviations. 
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OBJECTIVES 


The  Russell  (2006)  Ms  formula  has  opened  up  new  avenues  of  scientific  research,  such  as  the  development  of 
improved  regional  surfaee  wave  Q  models  (Stevens  et  al.,  2006;  Levshin  et  al.,  2006;  Cong  and  Mitehell,  2006)  that 
may  further  reduee  interstation  variance  of  the  magnitudes.  We  believe  that  application  of  the  A/v(VMAX)  technique 
to  Love  waves  is  the  next  logieal  step  in  the  scientific  process  that  could  lead  to  improved  discrimination.  Our 
objective  is  to  improve  seismic  event  screening  using  the  properties  of  Rayleigh  and  Love  waves.  We  are 
accomplishing  this  through  ( 1 )  the  development  of  a  Love-wave  magnitude  formula  that  is  complementary  to  the 
Russell  (2006)  formula  for  Rayleigh  waves  and  (2)  quantifying  differences  in  complexities  and  magnitude  variances 
for  earthquake  and  explosion-generated  surface  waves. 

RESEARCH  ACCOMPLISHED 


Love  and  Rayleigh  Wave  J/S(VMAX)  in  the  Middle  East 

We  have  applied  the  Ms  (VMAX)  analysis  (Bonner  et  al.,  2006)  using  both  Love  and  Rayleigh  waves  to  -100  events 
loeated  in  the  Middle  East.  A/v(VMAX)  is  for  estimated  both  Rayleigh  and  Love  waves  using  the  Russell  (2006) 
formula: 


Ms  =  log(aA)  +  ^log[sin(A)]+  0.003 1 


T  Vs 


A  -  log(/c )  -  0.43  -  0.66  log 


20 

T 


(1) 


Our  initial  hope  is  to  be  able  to  use  the  same  formula  for  both  phases.  The  details  of  the  processing  used  to  estimate 
A/v(VMAX)  are  described  in  Bonner  et  al.  (2006). 

The  study  area  (Figure  1 )  is  loeated  in  the  zone  of  continental  collision  between  the  Eurasian,  African  and  Arabian 
plates.  The  region  of  study  is  very  complex  and  spans  a  variety  of  different  tectonie  regimes.  The  seismieity  in  South 
and  Central  Iran  and  Turkey  is  in  the  upper  crust,  shallower  than  -20  km  (e.g.,  Engdahl  et  al.,  2006).  It  deepens 
toward  the  north  in  the  Alborz  region  in  Northern  Iran,  where  it  becomes  distributed  through  the  erust.  Further  to  the 
north  in  the  Central  Caspian  Sea,  the  seismieity  follows  the  Apsheron-Balkhan  Sill  and  reaches  depths  of  30-100 
km,  deepening  toward  the  north  Turkey  seismicity  is  dominated  by  strike-slip  focal  meehanisms  and  eoneentrated 
between  depths  of  10  and  20  km.  In  the  Zagros  fold-and-thrust  region,  most  of  the  earthquakes  are  shallower  than  30 
km,  w  ith  median  depths  of  15±7  km. 
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Figure  1.  Map  of  the  Middle  Eastern  events  for  which  MJVMAX)  for  Love  and  Rayleigh  was  estimated. 
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Wc  computed  A^VMAX)  for  over  100  seismie  events  loeated  in  this  region  with  reported  body  wave  magnitudes 
(mh)  between  3.8  and  5.6.  The  majority  of  the  location  and  magnitude  information  (w  ith  a  few  exceptions)  was 
obtained  from  the  NEIC  bulletin.  The  comparison  between  A/v(VMAX)  for  Love  and  Rayleigh  waves  is  shown  in 
Figure  2.  The  A^/v(VMAX)  computed  using  Love  waves  is  greater  than  the  magnitude  for  Rayleigh  waves  for  the 
majority  of  the  events  of  larger  magnitudes  (above  4).  For  smaller  events,  however,  we  observe  a  large  number 
of  events  with  the  Rayleigh  A/v(VMAX)  exceeding  the  Love  A/v(VMAX).  This  peculiarity  eould  be  eaused  by  either 
reduced  SNR  for  smaller  magnitude  events,  or  by  some  unknown  source  processes,  such  as  a  dominant  normal  fault 
mechanism.  In  addition,  regional  differences  in  the  wave  attenuation  and/or  anisotropy  eould  eause  changes  in  the 
amplitudes  for  the  rays  traveling  in  different  directions.  Sinee  the  station  coverage  is  not  homogeneous,  these 
propagation  effects  could  potentially  result  in  biases  in  A^v(VMAX)  estimate  for  smaller  events  with  limited 
sampling. 
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Figure  2.  Comparison  of  T/v(VMAX)  computed  using  Rayleigh  and  Love  waves.  Small  blue  circles  show 

earthquakes,  the  red  symbols  show  Ms  estimate  using  one  station  (BRVK)  for  some  Soviet  nuclear 
explosions. 

Azimuthal  dependency  of  the  estimated  Afv(VMAX)  for  several  larger  events  of  known  magnitudes  (mh)  is  shown  in 
Figure  3.  The  azimuthal  differences  are  most  likely  eaused  by  source  effects  (e.g.,  non-isotropie  souree  radiation), 
however  the  propagation  effects  (c.g.,  anomalies  in  regional  attenuation,  anisotropy)  can  play  a  certain  role.  The 
stations  located  between  the  azimuths  0°  and  60°  often  have  higher  magnitude  measurements  than  the  stations 
loeated  in  other  directions.  Notiee  that  for  larger  events  the  Love  M(VMAX)  is  greater  than  the  Rayleigh 
A/v(VMAX)  for  most  azimuths.  This  difference  is  more  pronounced  for  strike-slip  type  events,  and  beeomes  smaller 
for  the  thrust  events,  eommon  for  the  Zagros  region 

To  test  whether  the  azimuthal  differences  in  the  magnitudes  are  caused  solely  by  the  focal  meehanisms,  or  if  there 
are  propagation  effects,  wc  plotted  A/v(VMAX)  for  the  individual  stations  grouped  by  the  direction  of  propagation. 
Figure  4a  shows  the  cross-plot  between  Rayleigh  and  Love  A/v(VMAX)  for  the  stations  located  to  the  NE  from  the 
corresponding  events  (baek  azimuth  range  between  20°  and  70°).  There  is  a  significant  number  of  measurements 
with  Rayleigh  A/,(VMAX)  exceeding  Love  Afv(VMAX)  for  smaller  events.  The  measurements  for  the  stations 
located  to  the  NW  of  the  events  (azimuth  range  between  270°  and  360°)  show  smeared  distribution  for  broader  range 
of  magnitudes  (Figure  4b).  No  significant  anomaly  for  smaller  events  is  observed.  Therefore  it  is  likely  that  the 
reverse  in  Love  and  Rayleigh  A/v(V MAX)  for  smaller  magnitudes  is  caused  by  the  regional  propagation  effects, 
rather  than  by  the  source  eharaeteristies. 
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Figure  3.  Azimuthal  dependence  of  A/S(VMAX)  for  some  events  with  known  focal  mechanisms.  Blue  and 
green  circles  show  Rayleigh  and  Love  3/s(VMAX)  estimates,  respectively. 

Love  and  Rayleigh  Wave  A/S(VMAX)  in  the  Korean  Peninsula 

We  applied  the  A/v(VMAX)  technique  to  31  earthquakes  (Figure  5)  occurring  between  1996-2008  located  in  the 
Korean  Peninsula  and  surrounding  regions.  The  events  ranged  in  size  between  3.2  <  Mw<  5.1.  The  distances  to  the 
3C  stations  (red  circles  in  Figure  5)  recording  the  events  ranged  from  55  km  to  1900  km.  For  the  estimation  of  Love- 
wave  magnitudes  data  were  converted  to  transverse  motion,  while  the  Rayleigh-wave  magnitudes  were  estimated 
using  vertical  component  data.  The  analysis  resulted  in  298  single-station  estimates  of  A/v(Love)  and  266  estimates 
of  Ms(  Rayleigh). 
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Figure  4.  Comparison  of  A/S(VMAX)  computed  using  Rayleigh  and  Love  waves  for  individual  stations:  a) 
stations  located  to  the  NE  from  the  events  (back  azimuth  ranges  between  20°  and  70°);  b)  stations 
located  to  the  N\V  (back-azimuth  ranges  between  270°  and  360°). 
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Figure  5.  Map  of  the  Korean  Peninsula  region  events  (y  ellow  stars)  for  which  A/S(VMAX)  using  Love  and 
Rayleigh  waves  was  estimated.  The  red  circles  are  stations  used  in  the  analysis. 

Figure  6  provides  a  comparison  of  Love  vs.  Rayleigh  wave  A/v(VMAX)  estimates.  Applying  Equation  1  calibrated 
for  Rayleigh-waves  results  in  a  positive  bias  for  the  Love-wave  magnitudes.  This  result  is  consistent  with  the  results 
for  events  with  A/,(VM  AX)  >  4  in  the  Middle  East  study.  Only  6  of  the  3 1  network  magnitudes  had  a  larger 
Rayleigh-wave  magnitude  than  Love-wave  magnitude,  although  when  the  standard  deviation  was  considered,  these 
events  could  fall  below  the  line  representing  equivalent  Love-  and  Rayleigh-wave  magnitudes. 

In  addition  to  the  earthquake  dataset  we  analyzed  the  two  North  Korean  underground  nuclear  explosions  (UNE). 
Bonner  et  al.  (2008)  found  A/v(VMAX)=2.93  for  the  Rayleigh  waves  from  the  9  October  2006  North  Korean  UNE. 
There  were  no  Love  waves  registered  on  any  of  the  analyzed  stations.  The  25  May  2009  announced  UNE  in  North 
Korean  had  a  larger  magnitude  [A/S(VMAX)=3.6]  and  had  Love  waves  large  enough  for  analysis.  The  A/v(VMAX) 
estimated  for  the  Love  waves  was  3. 1 ,  which  placed  it  below  the  earthquake  population  in  Figure  6.  The  25  May 
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2009  event's  standard  surfaee  wave  magnitude,  based  on  Rayleigh  waves,  plots  above  the  Murphy  et  al.  (1997) 
screening  line  for  Mx:mb  and  is  perhaps  even  more  anomalous  than  its  predecessor  (Bonner  et  al.,  2006).  It  appears 
that  the  reason  the  first  event  did  not  separate  from  the  earthquake  population  well  was  not  due  to  eonvergenee  of 
the  populations  near  m y,  -4. 
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Figure  6.  Comparison  of  T/*(YfMAX)  computed  using  Rayleigh  and  Love  waves  for  the  Korean  Peninsula 
region.  Small  circles  show  earthquakes  while  the  star  shows  the  Ms  estimate  for  the  25  May  2009 
announced  nuclear  test. 
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Figure  7.  Comparison  of  network  A/S(Y\M AX)  estimates  for  Rayleigh-  and  Love-waves  versus  International 
Data  Centre  (IDC)  body  wave  magnitude  mh  for  Korean  Peninsula  earthquakes  and  the  two 
announced  North  Korean  nuclear  explosions.  The  dashed  line  represents  the  event  screening  line 
proposed  by  Murphy  et  al.  (1997). 
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Unequal  Variance  Assumptions  and  Improved  Discrimination 

We  have  attempted  to  use  the  potential  of  observed  variances  that  differ  significantly  in  earthquake  and  explosion 
A/v(VMAX)  populations  (Figure  7)  for  improved  event  screening.  We  have  considered  two  possible  methods  for 
incorporating  unequal  variances  into  the  discrimination  problem  and  compare  the  performance  of  various 
approaches  on  a  population  of  73  WUS  earthquakes  and  131  WUS  explosions. 

The  conventional  statistical  approach  to  separating  the  distributions  with  unequal  variances  would  be  to  derive  the 
likelihood  ratio  criterion  for  that  ease  whieh  yields  a  hybrid  discriminant  composed  of  a  linear  term  and  a  quadratic 
term  with  weights  dependent  on  the  unequal  eovarianee  matrices  for  the  two  populations.  Taking  the  solution 
involving  the  likelihood  ratio  guarantees  that  the  probability  of  eorreetly  ruling  out  an  event  as  an  explosion  subjeet 
to  a  fixed  probability  of  failing  w  ill  be  maximized  under  the  multivariate  normal  assumption  for  the  input  variables, 
usually  taken  to  be  Ms  and  w/,.  Sometimes,  additional  tuning  parameters  are  added  to  eovarianee  matrices  that  may 
improve  performance  in  particular  samples  or  w  hen  the  joint  Gaussian  assumptions  are  not  satisfied  (see  Friedman, 
1989;  Anderson  et  al.,  2007).  The  eommon  assumption  that  the  eovarianee  matrices  of  the  two  populations  are 
approximately  equal  reduees  the  classification  function  to  a  linear  function  of  the  input  variables.  For  example,  the 
simple  difference  M-mh  is  often  used  as  a  further  approximation.  This  approach,  however,  is  not  recommended 
beeause  of  poor  performance  in  sample  populations  sueh  as  the  one  considered  here. 

An  alternate  method  proposes  replacing  the  Ms  component  by  Ms  +  a*a,  where  a  denotes  the  interstation  standard 
deviation  obtained  from  the  stations  in  the  sample  that  produced  the  Mx  value.  In  this  eontext,  we  interpret  this  as 
replacing  the  usual  linear  discriminant  a*A/v  +b *mh  by  a*A/v  +b*w/,+e*o  (Figure  9).  We  also  estimate  the  optimum 
hybrid  linear-quadratie  discriminant  function  resulting  from  the  unequal  variance  assumption.  While  the  input 
standard  deviations  will  not  be  normally  distributed  (they  follow  the  ehi  distribution),  the  linear  approximation  may 
be  reasonable  for  the  same  reasons  that  the  linear  discriminant  function  works  in  the  usual  case.  While  the  two 
discriminant  functions  resulting  from  the  two  theoretical  interpretations  of  the  unequal  variance  function  did  slightly 
better  in  the  test  samples  used  here,  all  methods  exeept  Ms  - m /,  did  extremely  well  using  the  WUS  dataset.  We  plan 
on  evaluating  this  method  further  using  our  entire  dataset  of  earthquakes  and  explosions. 


Figure  8.  Frequency  histograms  for  the  station  standard  deviations  (top)  and  variances  (bottom)  for  73  WUS 
earthquakes  and  131  explosions. 
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Figure  9.  Conventional  multivariate  discriminant  analysis  combining  Ms(VMAX)-wa  (left),  MS(V MAX),  mb 
and  standard  deviation  (std)  (middle),  and  multivariate  linear-quadratic  discriminant  analysis 
(right)  for  73  WUS  earthquakes  and  131  explosions. 

A  New  Metric  for  Intrastation  Complexity 

We  have  also  studied  the  complexity  of  the  “magnitude  spectra”  at  eaeh  station.  Our  hypothesis  is  that  explosion 
speetra  should  have  fewer  foeal  meehanism-produeed  complexities  in  the  magnitude  speetra  than  earthquakes 
(Figure  10).  We  have  tried  several  different  methods,  including  the  intrastation  std  as  a  function  of  period  shown  in 
Figure  10  and  a  “differencing”  approach.  For  periods  /  between  0  and  25  seeonds,  we  estimate: 

AXfs  =  Ms(i)-Ms(i+\).  (2) 

We  start  at  i-9  seeonds  to  minimize  edge  effeets.  This  results  in  n=\7  new  magnitude  differences,  zfWv,  for  eaeh 
station.  It  may  be  fewer  if  estimates  at  some  periods  fail  a  signal-to-noise  ratio  test. 

We  then  estimated  several  different  metries,  M,  from  4v/s  for  comparison  of  the  differenced  magnitudes.  Thus  far, 
we  have  settled  on: 


M  = 


n 


(3) 


Figure  1 1  shows  the  histograms  for  the  earthquake  and  explosion  populations  M while  Figure  12  shows  the 
cumulative  distribution  functions  of  M  The  metrie  M  works  well  for  the  application  of  trying  to  separate  complexity 
of  souree  speetra.  At  90%  explosion  eonfidenee,  around  .02  M r,  there  is  about  a  30%  earthquake  eonfidenee. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Preliminary  research  shows  that  Love-wave  magnitudes  for  earthquakes  are  often  equal  to  or  larger  than 
Rayleigh-wave  magnitudes.  Conversely,  for  explosions,  Love-wave  magnitudes  are  typically  smaller  than  Rayleigh 
wave  estimates,  or  below  background  noise  levels.  However,  we  observe  a  number  of  smaller  events  with  reversed 
pattern.  Interestingly  enough,  this  peculiarity  is  observed  for  the  Middle  East  dataset,  but  not  for  Korean  dataset.  We 
will  continue  to  examine  this  phenomenon,  possibly  incorporating  Ms  maximum  likelihood  estimates,  with  hopes  of 
porting  the  results,  theory,  and  statistical  p-values  into  the  Event  Classification  Matrix  (ECM;  Anderson  et  al., 

2007).  Additional  aspeets  of  surfaee  wave  propagation,  ineluding  a  new  intrastation  complexity  metrie  and 
differences  in  explosion  and  earthquake  magnitude  variances,  also  show  promise  for  improved  event  screening. 
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We  will  also  develop  an  improved  formula  for  estimation  of  Aft(VMAX)  using  Love  waves.  Equation  1  was 
developed  using  empirical  relationships  using  the  Rayleigh  waves,  which  may  or  may  not  be  calibrated  for  Love 
waves. 
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Figure  10.  A  proposed  ne>v  approach  for  improved  discrimination.  (Left  Column)  Traditional  A/S(V MAX) 
w  here  the  circle  is  the  period  of  maximum  amplitude.  (Right  Column)  The  mean  of  a  station's 
estimate  and  its  standard  deviation  are  formed.  Our  results  suggest  that  the  intrastation  standard 
deviation  is  typically  less  for  explosions  than  earthquakes  for  our  VVUS  database. 
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Figure  1 1.  Comparison  of  the  new  intrastation  metric  Aw,  for  NTS  explosions  (left)  and  VVUS  earthquakes 
(right). 
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Cumdrtwt  DtttnbiXion 


Figure  12.  Cumulative  distribution  functions  for  metric  M  (Equation  3;  Figure  1 1)  for  VVUS  earthquakes  and 
NTS  explosions. 
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